On Scaling Solutions with a Dissipative Fluid 
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We study the asymptotic behaviour of scaling solutions with a dissipative fluid and we show that, 
contrary to recent claims, the existence of stable accelerating attractor solution which solves the 
'energy' coincidence problem depends crucially on the chosen equations of state for the thermody- 
namical variables. We discuss two types of equations of state, one which contradicts this claim, and 
one which supports it. 

<n : 

Current observational evidence suggests that the Universe is fiat and accelerating (T|. Since the average energy- 
density of matter is relatively low there should consequentially exist a dark energy smoothly distributed with roughly 
the same magnitude as the matter energy. In general this requires fine-tuning of the initial conditions and has become 
known as the 'energy' coincidence problem. The simplest explanation for this dark energy is a vacuum energy density; 
however, due to the discrepancy between the "observed" value of the cosmological constant and that predicted by 
particle physics alternative mechanisms should be considered. One model that has recently attracted attention 
explains the "missing" energy in terms of the existence of a minimally coupled scalar field ("quintessence" field) [3|. By 
a suitable choice of the scalar field interaction potential, the scalar field evolves becoming the dominant contribution 
to the energy density today. These type of solutions (matter scaling solutions) were first studied in Q and their 
t— I . stability, for the exponential potential case, were analyzed in ||. However, since FLRW solutions with a perfect fluid 
and a scalar field evolve towards a scaling solution that does not accelerate (unless matter violates the strong energy 
condition), Chimento et al. Q considered a dissipative fluid by means of the existence of a bulk dissipative pressure. 
Bulk viscosity may be associated with particle production Q , may be understood as frictional effects that ap Dears in 
£^ mixtures or even to model other kinds of sources (e.g. the string dominated universe described by Turok [10| or a 

■ scalar field p|). 

{NJ | In H the authors studied the effect of the bulk viscosity by using the truncated Israel-Stewart theory Q and claimed 
that there exists a stable attractor that has an accelerated expansion thereby solving the energy coincidence problem. 
""^ | However, their analysis was not complete: the differential equation for the function gamma defined in (10) of their 
q-i article was not taken into account and the equation of state for the bulk pressure was not stated. In this paper we 
study the asymptotic behaviour of a FLRW model with a dissipative fluid and a scalar field with exponential potential 
using two different equations of state for the bulk pressure. 

We will show that the claim stated in || , that there exist stable accelerating matter scaling solutions with non-zero 
bulk viscosity, cannot be substantiated. We will show that the conclusions in |6) depend crucially on the chosen 

■ equation of state for the coefficient of bulk viscosity, £, and the relaxation time, t. 
H ' We assume a FLRW metric: 



in 



C3 



ds 2 = -dt 2 + a 2 {t) dn 2 (1) 
with a source given by a fluid with bulk viscosity and a scalar field with a exponential potential: 

T ab = T* b +T™ b , (2) 

where 



T* b = <f>a<t>b-9ab fyc^ + V^ 



(3) 



and 

T ab = pu a u b + {p + <7){u a u b + g ab ), (4) 
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where a is the viscous pressure, p is the thermodynamic pressure and p is the energy density. Einstein's equations 
arc then: 



H 



-H 2 - hp + 3p + 3a + 2<j> 2 - 2V) 





3H 2 = p+ l -tf 



n 



p = -3H(p + p + a) 

where H = a/a. 

The viscous pressure a is assumed to satisfy the truncated Israel-Stewart equation Q: 

a + ra = -3£H, 



(5) 

(6) 
(7) 

(8) 



(9) 



where £ is the coefficient of bulk viscosity and r is the relaxation time (£ > 0, r > 0). Strictly speaking, eq. (9) is valid 
only when the fluid is close to equilibrium; however, we assume its validity even when the fluid is far from equilibrium. 
(A more rigourous and comprehensive analysis would utilise the full transport equation of the Israel-Stewart theory.) 
We will assume a linear barotropic equation of state p = (7 — l)p, where the constant 7 satisfies < 7 < 2 and two 
different relations for the coefficient of bulk viscosity and relaxation times. 
In the first case we consider the following relations: 



which were introduced by Belinskii et al. ]T^ |. 

We now introduce a new set of normalized variables: 



p 



ap 



(10) 



n = 
r = 



3H 2 ' 
1 V 



and a new time r defined by: 



h = H x - 2m , 
H(t) dt = dr. 



= _ 

V6H' 

m ^ 77' 



The Friedmann equation (6) now becomes: 



l-17 2 -* 2 -r = - 



K 



H 2 a 2 ' 



and Einstein's equations are now: 



Q' = n(-3-y + 2x) - E 
£' = -9fi-F 



(3fi)( 1 " m )/i + 2a 



*' = *(x-3)-^|r 



— (1 — 2m)hx 
where ' denotes derivative with respect to r and 



r' 

ti 



(11) 

(12) 

(13) 

(14) 
(15) 

(16) 

(17) 
(18) 



-(37-2)fi- 



2* 2 -r. 



(19) 



When m = 1/2 the equation for H decouples from the other equations and the system becomes a four-dimensional 
system. We will not consider this particular bifurcation situation just yet as it is a special case of the next equation 
of state considered below. 
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Before studying the equilibrium points of the above system we present some relevant expressions. The energy 
density of the scalar field is given by: 



P4> = \^ 2 + V = 3H 2 (V 2 + T) 



and 



The adiabatic index of the scalar field is: 



'AH 2 



+ r. 



70 



P<p+P<f> 
P<t> 



2* z 



U 2 + v * 2 + r' 



The deceleration parameter is given by: 



q = x — 1. 



(20) 



(21) 



(22) 



(23) 



We are looking for those equilibrium points that have non vanishing matter and scalar field; i.e., ft ^ 0, \I> ^ and 
f^O. From equation (17) one such point corresponds to: 



x = -k — 
2 



(24) 



Equation (18) implies that h = and therefore H — > when to < 1/2 and H — > oo when m > 1/2, as the point is 
approached. From equation (16) we obtain: 



k 



From equation (15) we obtain: 

= -9Q + 2.tS, 

Substituting this expression in equation (14) we get: 

V6 



9 



(25) 



(26) 



(27) 



Since T > 0, equation (25) implies that \& must be negative and therefore only the solution with the minus sign is 
valid in the above expression. Furthermore, it is easy to see that the condition '5 < implies that 7 < 3/2. Finally, 
to obtain £ we substitute equations (24)-(27) into equation (19) and we obtain: 



£ = - 



9 1 
V6fc* 



(28) 



So the equilibrium point is given by: 



n = 1 

E 



r 

ft = 



7 + vV + 4 

(7+^Tl) 

3 
4fc2 

3 



l-^(7+V^T4) 



(2 - 7) (7 + vV + 4) - 2 



7< 2- 



(29) 
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The energy density of the scalar field is 

3 



n = ^ 2 + r = — ( 7 + V7 2 + 4) , n + n* = i, (30) 



and the adiabatic index is 

2* 2 



l<t> = 



5 (7 + vV + 4) ^ 7 - (31) 



* 2 + r 

The deceleration parameter in this case is given by: 

9 = ^(7+V / 7 2 +4)-l>0. (32) 

As in the matter scaling solutions, this equilibrium point represents a flat universe in which the energy density of 
the scalar field scales with that of the matter; however, in these solutions, the adiabatic indexes are different. More 
importantly: this solution is not accelerating. 

The stability of this equilibrium point is determined by the sign of the eigenvalues of the matrix of the linearized 
system around the point. After a lengthy calculation one can find that three of the eigenvalues are: 

Sv/t^+4, lil+V^FTl) -2, ^( m _I)( 7+v /^T4). (33) 

The sign of the two first eigenvalues are positive whereas the sign of the third depends on the value of m. Thus this 
equilibrium point is not stable. The signs of the two other eigenvalues are difficult to determine. 

It is important to note that this is the only equilibrium point with ^ 0, $ ^ and r ^ 0. If we relax the 
last condition we obtain, in addition to the former point, a set of equilibrium points corresponding to T — 0; i.e., a 
massless scalar field. From equation (16) we get x — 3 and from (18) h = 0. From equations (14) and (15) we obtain: 

£=^0, 7=| (34) 

And finally from equation (19) we obtain: 

f) + * 2 = l. (35) 

In this case 7^ = 2 and fl^ — \& 2 , (fi + fi^ = 1). Again this point is not stable since the eigenvalues of the matrix of 
the linearised system are: 

0, 6(m--), 6 + V6fc*, 4, 15. (36) 

We shall now consider the second case by using the following equations of state: 

£ = £op m H l - 2m , T = T p m - r H 2r - 2m - 1 , (37) 

as used in Q . 

The dynamical system is unchanged except for the evolution equation for the viscous pressure, which now becomes, 
in terms of the new variables ( [Tl"| ) 

£' = -aCf + S {—bn r ~ m + 2x] (38) 



where 



nl+rt or~m 

a = - 2°, and b=- (39) 

TO TO 

replace our original positive constants £0 an< 4 T o- As in ||l3f we take r = 1, and we take to < r to ensure continuity 
of d38| ) . Note that as H does not appear in this equation the equation for h decouples from the rest of the dynamical 
system and we need only study the resulting four-dimensional system. Note also that when to = 1/2, we have a slight 
generalisation of the equation of state used above. 
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There are two equilibrium points which satisfy the requirement 17 > 0,T ^ 0, and ^ ^ 0. The first has 

E = -17 (/c 2 (17 - 1) + 3 7 ) (40) 

* = -^fc(n-i) (4i) 

r = (q - 1) ( -k 2 (i - n) - i ) , (42) 



6 



where 17 6 (0, 1) must satisfy 
From @) and T > we see that 

hence, from (E3T), we see that 



(k 2 (17 - 1) + 3 7 ) (617 1 -™ + A: 2 (17 - 1)) - a = (43) 



fi>l + ^; (44) 



a < 3 (7 + 2) (6 + (45) 



must be satisfied for this point to exist. 

The scalar field has density parameter and adiabatic index 



The deceleration parameter is given by 



17 = 1 - 17, 70 = \k 2 (1 - 17) . (46) 



Q = \k 2 (1 - fi) - 1 = §7* - 1. (47) 



This equilibrium point will be stable if all the real parts of the eigenvalues of the Jacobian (of the rhs of the 
dynamical system) are negative. One of the eigenvalues A is easy to determine: 

A = fc 2 (l - 17) - 2 = 2g; (48) 

thus if this point is not inflationary it will not be stable cither. The remaining eigenvalues may be determined but 
the resulting expressions are extremely large. Therefore we will present the results for this point graphically, in 
figures (|l|) - (Q). As we only wish to qualitatively show the different behaviour which may occur, we only show a 
few representative cases, and show the variation of behaviour with the parameter a. (Recall that tq — 3 1 ~ m /b, and 
£o = a/3 1+m 6.) The effect of increasing b may be seen from the figures: for small b (Fig. [I]) the upper curve (high 
density solution) is stable (as a is increased) until 17 = 1 is reached at a ~ 3 [for a larger than this the high density 
solution no longer exists; cf. Eq. (f45|)]; similarly for b = 5 (Fig. ||), once the upper curve becomes stable (as a is 
increased) at a ~ 2, it remains so until Q = 1, which happens in this case at about a ~ 14. For the case b = 10 
(Fig. ||), the higher density solution becomes stable at a ~ 6, and fi = 1 is reached at a ~ 30. The lower density 
solution appears never to be stable. 

In the figures, we have arbitrarily set 7 = 1, m = 0.5, and k — 2. The effect of increasing 7 or m is to push the 
curves down to lower 17 (the lower density solution eventually drops off the plots); this does not change the conclusions. 
Lowering the value of m allows the solution to become stable for smaller a (while pushing the curves up to higher 17); 
e.g., in figure (||) with m = 0.1, the upper curve becomes red (stable) before a = 1. The effect of changing k is shown 
in figure (Q). 

The important point here is that stable solutions exist for a large non-trivial region of the five-dimensional parameter 
space {a, 6, to, k, 7} of the models. 

The other equilibrium point of interest here is characterised by the following (when to =/= 1) : 



k a + 2(37-2), 

E = -17(37- 2), (50) 

* - <-» 

r = 5F- < 52 > 
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FIG. 1: A plot of f2 versus a for the parameter values k — 2, m = 0.5, 7 = 1 and 6 = 1, obtained as the solutions of (|43|). The 
thin grey curve (lower) represent solutions which are unstable, while the thicker black curve represents the stable solution, and 
is thus inflationary, by (ffq). 

In 




For this point 

2 2 

7< ^ = 3' q = - ( 53 ) 

Therefore this point is not inflationary. However, this point in general has non-zero curvature. The stability of this 
point is difficult to determine analytically, so again we present the results graphically in figure (||) . Again we see that 
there exist stable solutions. 

This graph (||) suggests that, for these parameter values, all solutions with fl > 0.5 are unstable, a fact implied 
by further investigations (although the limiting value of depends on k: indeed for high k there are no unstable 
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FIG. 3: A plot of Q. versus a for the parameter values k = 2, m = 0.5, 7 = 1 and 6 = 1, 5, 10, obtained as the solutions of (^3J). 
This is an amalgamation of the previous two plots, with the additional b — 10 case represented in black, with the previous 
cases given in shades of gray, as shown in the key. Again, the thick curves represent the stable solution, and the thinner curves 
the unstable solutions. The case 6 = 10 only becomes stable at high a. 
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FIG. 4: A plot of f2 versus k for the parameter values a = 1, m = 0.6, 7 = 1 and 6 = 1, obtained as the solutions of (|l3[). 
The effect of increasing k alters the value of Q, but not the stability; the stable solution (the thick black curve) exists for all k, 
while the unstable solution (thin, grey) becomes available for k > 2. 



solutions). This then suggests that the stability of the solution depends on the curvature: for fl < ft,/, = 2/ k 2 (= 0.5 
in figure (||)) we have the stable solutions, and these have negative curvature, by ([l3|). 

For the previous two points we have demonstrated the existence of stable solutions; however, we ignored the 
traditional thermodynamical assumption that the magnitude of the viscosity must be less than the isotropic pressure 
(we neglected this in the figures by setting 7 = 1). However, we note that this does not affect the result. (So, e.g., in 
the latter equilibrium point this requirement would set i < 7 < 3; stable solutions exist in this case too.) 
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FIG. 5: A plot of f2 versus a for the parameter values k = 2, m = 0.6, 7 = 1 (so Q<£ = 0.5); three values for b are represented, 
6 = 5, 10, 15, with b = 5 being the upper curve, and 6 = 15 the lower. The thin grey curves represent solutions which are 
unstable, while the thicker black curves represent the stable solutions. 



Finally, if we relax our restriction that T ^ 0, we have another two equilibrium points, given by (when ra/1) 

n.USfc^Y*, (54 ) 



^+18(7-2), 

E = 30(2-7), (55) 
* 2 = l-n = fV (56) 

For this point q = 2, and 7^ = 2, so it is not inflationary. Three of the eigenvalues are 

4, 3(2m-l), 6 + V6fc*; (57) 

hence this point is not stable. 

We have shown that the equations of state play a crucial role in the final evolution of non-perfect fluids with a scalar 
field. Therefore the claim that a dissipative fluid minimally coupled with a scalar field can resolve the coincidence 
problem, although suggestive, depends dramatically of the equations of state and on the parameters associated with 
them. 
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